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A new apparatus for measuring the thermal conductivity of fluids is described. This is an absolute method 
utilizing a transient hot wire. Measurements are made with a 12.7 /tm diameter platinum wire at real times of up 
to 1 second. The data acquisition system includes a minicomputer and a digital voltmeter. The experimental 
core of the system incorporates a compensating hot wire in a Wheatstone bridge circuit. The cell containing the 
core of the apparatus is designed to accommodate pressures from to 70 MPa and temperatures from 70 to 320 
K. Oxygen was measured over a wide range of physical states including the dilute gas, the moderately dense gas, 
the near critical region, the compressed liquid states, and the vapor at temperatures below the critical temper- 
ature. Performance checks of the apparatus were conducted with nitrogen, helium and argon. Measurement of 
rare gases allows a direct comparison to the kinetic theory of gases through the viscosity. A second check looks 
at the variation of the measured thermal conductivity as a function of the applied power. The precision (2 a) of 
the new system is between 0.5 and 0.8 percent for wire temperature transients of 4 to 5 K, while the accuracy is 
estimated at around 1.5 percent. 
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1. Introduction 



2. Method 



Thermal conductivity values are necessary whenever a 
heat transfer problem is to be evaluated. In addition, ther- 
mal conductivity is a property of fundamental interest in 
developing the theory of fluids. Accurate measurements of 
thermal conductivity are of considerable difficulty. Methods 
and geometries abound each with its adherents and its in- 
herent drawbacks. The steady state hot wire experiment is 
one of the older, well established methods. The transient 
hot wire method, however, has come into its own only with 
recent advances in digital electronics. The evolution of the 
modern transient hot wire experiment can be traced from 
the early experiments of Pittman [l] 1 , Haarman [2], and 
Mani [3]. However, the major exposition of both theory and 
application of this method was accomplished by Kestin and 
his coworkers during the last decade [4-9] for the gas phase. 
Similar instruments have been used by Wakeham and his 
coworkers for both gas and liquid measurements [10, 11] 
and by de Castro and his coworkers for liquid meas- 
urements [12-13]. 
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A hot wire system normally involves a vertical, cylindrical 
symmetry where the wire serves both as heating element 
and as thermometer. Almost without exception platinum is 
the wire of choice. The mathematical model that one at- 
tempts to approximate is that of an infinite line source of 
heat suspended vertically in an infinite medium. The 
method is labeled transient because the power is applied 
abruptly and the measurement is of short duration. The 
working equation is based on a specific solution of Fourier's 
law and can be found in standard texts (see for example 
reference [14], page 261). 



m -T r . f = AT = ^- ln ^ 







(1) 



Where T[t) is the temperature of the wire at time t; 

T r .f is the reference temperature, the temperature of the 
cell; 

is the applied power; 

is the thermal conductivity of the fluid, a function of 
both temperature and density; 
is the thermal diffusivity of the fluid, i.e., K = 
\l qC p . K is normally taken at the temperature T r . f 
and is nearly constant since the fluid properties do 
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not change drastically with a small increase in 
temperature; 

a is the radius of the wire; and 

in C= 7 where 7 is Euler's constant. 
The relation given by eq (1) implies a straight line for a plot 
of AT versus (n(t). In practice systematic deviations occur at 
both short and long times. However, for each experimental 
measurement there exists a range of times over which eq (1) 
is valid, that is the relation between AT and tn{t) is linear. 
This range of validity is determined from 250 measured 
AT-t pairs by selecting a beginning time h and an ending 
time t 2 . The slope of the A T versus ln(t) relation is obtained 
over the valid range, i.e., between times fj and t 2 , and using 
the applied power the thermal conductivity is calculated 
from eq (1). The temperature assigned to the measurement 
of X is given by 

T= T,., + ^[ATitJ + A1\t 2 )} (2) 

The density assigned to the measurement of X is taken from 
an equation of state using an experimentally measured 
pressure and the temperature assigned above. The equa- 
tions of state used for nitrogen and helium are given in 
subsequent sections. The experimentally determined tem- 
perature rise of the wire is AT W . A number of corrections ac- 
count for the departure of the real instrument from the 
ideal model: 



AT = AT W - E 5T, 



(3) 



These corrections 8T { have been fully described elsewhere 
[5]; the most important at lower times is 57\, the effect of 
the finite heat capacity of the wire. 

3. Apparatus 

The various elements of the apparatus, the hot wires, the 
high pressure cell and wire supports, the Wheatstone 
bridge, the cryostat, the measuring and control circuitry, 
the sample handling system, and the minicomputer are de- 
scribed here and in figures 1-6. 

3.1. The Hot Wires 

Platinum wire is the normal choice because the resist- 
ance-temperature relation of platinum is well known [15]. In 
the resistance thermometer grade of the wire the smallest 
sizes available commercially are 7 /im, 12.7 pm, and 25.4 /*m 
in diameter. The smallest of these was considered to be too 
fragile for the present application, the largest would not 
provide a sufficiently large resistance for the cell size under 
consideration. Therefore, the 12.7 ftm diameter wire was 
chosen for both primary and compensating hot wire. The 



long or primary hot wire is approximately 10 cm in length. 
Its resistance varies from about 20 Q at 76 K to 90 Q at 298 
K. The short or compensating wire is approximately 5 cm in 
length and its resistance varies from 10 to 45 Q. To allow for 
thermal expansion over the range of operating tempera- 
tures the wires are mounted with a small weight attached at 
the bottom. The wires are mounted in the hard drawn (as 
received) state, because in the annealed state platinum has 
a rather low tensile strength [16]. 

3.2. High Pressure Cell and Wire Supports 

The cell has been designed to withstand pressures of up 
to 70 MPa. In addition^ if the cell is to be isothermal, its 
thermal conductivity should be high. Accordingly, the cell is 
made from beryllium copper, a material which meets these 
contradictory requirements best. Cell sizing, cell closure, 
and wire support are shown in figure 1. The cell is quite 
similar to those used in other experiments in this laboratory 
[17-20]. It is mounted with the closure facing down. Three 
solid steel pins serve as electrical leads. They pass through 
the high pressure closure (stainless steel) and are insulated 
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Figure 1. The high pressure cell and wire support unit. 
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from it using polymide ferrules [21]. The capillary which ad- 
mits the sample is soldered into the top of the cell. Much of 
the interior of the cell is filled by the wire support unit 
which is made from copper and includes a centerpiece, two 
covers, a top and a bottom cap. The wire support unit is 
mounted on the high pressure closure. When assembled the 
wire support unit provides two cylindrical wells 9 mm in 
diameter and 11 cm long to accomodate the hot wires. 
Mounted in the top cap with friction fit are two 5 mm sec- 
tions of a small coaxial cable. The hot wires are soft 
soldered to the center leads of these coaxial sections. To the 
bottom of the long hot wire is soldered a short section of 
teflon covered wire which terminates in a 4-80 brass nut, 
the weight. The teflon covered section provides isulation be- 
tween bottom cap and the extension of the long hot wire as 
well as centering for it. The long hot wire section is con- 
nected to one of the steel pins by a small loop of 12.7 /an 
copper wire. The short hot wire is soft soldered to the other 
coaxial section in the top cap. To its bottom end a short 
length of copper wire centered in a second 4-80 brass nut is 
soldered. The leads on the short hot wire side are completed 
by a loop of 12.7 /tm copper wire which connects to a 5 cm 
length of coaxial cable. This section of coax is friction fit 
into the bottom cap and is sufficiently rigid to stand alone. 
A second loop of copper connects this coax to the second 
steel pin. Long and short hot wire sections are connected 
together above the top cap with the center tap soldered in 
the middle of this connection. The center tap is a 11 cm 
long section of the coax which at the bottom end is con- 
nected to the third steel pin. 

Liquid oxygen safety is one of the additional design con- 
siderations for the cell since the interior of the cell will be 
exposed to very high pressure 70 MPa (10,000 psi) liquid. 
The materials directly exposed to liquid oxygen have been 
limited to beryllium copper, copper, stainless steel, silver, 
teflon, and a polyimide (kapton) all of which have been 
found to be "oxygen compatible" [22]. Cleaning procedures 
for cell, wire supports, capillary and sample handling 
system were extensive [23]. 



measured directly as a function of time with a fast response 
digital voltmeter (DVM). The DVM is controlled by a mini- 
computer which also handles the switching of the power and 
the logging of the data. The automation of the voltage 
measurement follows the work of Mani [3] who used a 
similar arrangement with a transient hot wire cell to meas- 
ure resistance by the four lead technique rather than using 
a bridge. 

Figure 2 shows the Wheatstone bridge circuit. Each arm 
of the bridge is designed to be 100 Q, two arms /?j and R 2 
are standard resistors. The resistance in each of the other 
arms R 3 and R t is a composite of the hot wire, the leads into 
the cryostat and an adjustable ballast resistor. The leads are 
roughly 6 Q at room temperature and 2 Q when the cell is at 
76 K. The ballast resistors allow each working arm to be ad- 
justed to a value of 100 Q. 

The measurement of thermal conductivity for a single 
point is accomplished in two phases. In the first phase the 
bridge is balanced as close to null as is practical. To start 
with, switch 1 is turned from dummy to the bridge while 
switch 2 is open. With a very small applied voltage, 0.1 
Volts normally, and the cell essentially at constant tempera- 
ture, the voltages are read on channels through 7. The 
lead, hot wire, and ballast resistances are calculated from 
the ratios of the appropriate channel voltgage to the voltage 
across the standard 100 Q resistor on channel zero. The 
ballasts are adjusted until each leg is approximately 100 ft. 
Finally, with switch 2 closed, the bridge null is checked on 
channel 6. The second phase incorporates the actual ther- 
mal conductivity measurement. The power supply is set to 
the applied power desired, switch 2 is closed, and switch 1 is 
switched from dummy to bridge. The voltage developed 
across the bridge as a function of time is read on channel 6 
and stored. The basic data is a set of 250 readings taken at 3 
ms intervals. Finally the voltage on channel is read to de- 
termine the exact applied power, and the power is switched 
back to the dummy resistor. 

3.4. The Cryostat 



3.3. The Wheatstone Bridge 

Precision measurements of resistance can be made by us- 
ing a four lead technique or by using a Wheatstone bridge. 
In the present apparatus we follow the general development 
of the hot wire instrument pioneered by Haarman [2], de 
Groot, et al. [4], Assael, et al. [10] and de Castro, et al. [12] 
and use a Wheatstone bridge to measure resistances. End 
effect compensation is provided by placing the long hot 
wire in one working arm of the bridge and a shorter, com- 
pensating wire on the other. In contrast to other instru- 
ments where values of time are measured in the present 
instrument the voltage developed across the bridge is 



The cryostat for the apparatus, shown in figure 3, is 
adapted from a general design first used in a PVT appa- 
ratus at this laboratory [24]. The cell is connected by the 
reflux tube to the inner refrigerant tank. The reflux tube is 
filled with gas which sometimes corresponds to the liquid 
used as the refrigerant. Varying the gas pressure changes 
the amount of refrigeration applied to the cell. A cooling 
yoke and ring insure that the refrigeration is applied 
primarily at the bottom of the cell. In this way a slight gra- 
dient can be maintained between cell top and cell bottom. 
To avoid convection inside the cell we normally maintain 
the cell bottom slightly colder than the cell top. The cell and 
its radiation shield is located in a vacuum environment. A 
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FIGURE 2. The Wheatstone bridge circuit. 



tube which serves as the inner vacuum line also admits lead 
wires and the capillary into the inner vacuum space sur- 
rounding the cell. The leads are tempered on two rings, one 
is attached directly to the inner refrigeration tank, the other 
is the guard ring which is attached to the reflux tube. The 
cell radiation shield in turn is attached to the guard ring. 
The inner refrigerant tank is insulated by the outer vacuum 
system, which in turn is protected by the outer refrigerant, 
usually liquid nitrogen, in a dewar. All of the various access 
tubes, fill and vent lines, and vacuum lines pass through the 
mounting plate which rests on the mounting bracket at- 
tached to the wall. Cell alinement is achieved by adjustment 
of three set screws. 

3.5. Cryostat Temperature Control 

A diagram of this circuitry is shown in figure 4. The cell 
temperature is monitored with the platinum resistance ther- 
mometer which is mounted in the cell. The PRT resistance 
is measured by the four lead method as a ratio against a 10 
standard using a microvolt potentiometer, a stable cur- 



rent supply, and a high gain dc null detector. The amplified 
output of the detector is used in a feedback loop to provide 
power to the cell heater [25]. The ratio of power applied to 
cell top and cell bottom is monitored with a thermocouple 
and adjusted by hand. The temperature of the guard ring 
and cell radiation shield are maintained close to the cell 
temperature by the use of a second feedback loop. This loop 
includes a differential thermocouple and a low level dc 
voltage detector whose amplified output is routed to a 
power supply which in turn feeds the ring and shield heater. 
A separate, manual heater is mounted on the capillary. This 
heater is operated intermittently as needed. Its function is 
to ensure that the sample does not freeze in the capillary 
should by chance some unusual temperature conditions pre- 
vail at some spot in the capillary. 

3.6. Sample Handling and Vacuum Systems 

Block diagrams for both systems are shown in figure 5. 
The vacuum system is conventional except for an automatic 
back pressure control with alarm and solenoid valve. These 
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items were installed as a safety precaution should the silver 
plated copper O-ring in the cell closure spring a leak and 
discharge oxygen into the vacuum space. 

The sample handling system was also designed with oxy- 
gen in mind. It consists of the sample from a gas bottle, a 
molecular sieve to capture moisture, and a small diaphragm 
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FIGURE 3. The cryostat. 



compressor as pressure intensifier. The compressor uses an 
oxygen compatible oil and has specially hardened check 
valves and seals. A pressure relief valve was provided in the 
cell pump-out line to protect against the possibility of an ac- 
cidental oxygen overpressure, even though the forepump is 
charged with an oxygen compatible pump oil. With a few 
changes the sample system was able to handle a liquid sam- 
ple, propane. Ordinarily the compressor requires an input 
pressure of about 5 MPa (700 psi). However, if liquid is pro- 
vided at the intake and sufficient time between strokes is 
allowed, then the compressor handles a liquid sample quite 
well. For propane, removing the molecular sieve and turn- 
ing the supply bottle upside down ensured a direct flow of 
liquid. 

3.7. The Minicomputer 

A simplified block diagram of the minicomputer and its 
peripherals is shown in figure 6. The CPU uses 18 bits per 
word with BASIC as the operational language. Program 
storage is on floppy disks, program input on the CRT, pro- 
gram listings are handled on the printer. Data can be dis- 
played on the CRT and on the printer, and can be stored on 
floppy disk or magnetic tape. Input voltages, i.e., voltages to 
be read are routed from the 50 input channels by the multi- 
plexer to the DVM where the analog to digital (A/D) conver- 
sion is effected. Output voltages, i.e., voltages to be used for 
experimental control are processed by the logic control unit 
where the digital to analog conversion takes place. These 
voltages are available at any of six output channels. 

4. Data Measured 

In the course of making a single thermal conductivity 
measurement a large number of variables are measured and 
recorded. The minicomputer program which controls the 
measurement is shown in appendix I. An example of a data 
file as assembled in the minicomputer and then transferred 
to magnetic tape is shown in appendix II. The first two lines 
of the data file are keyed in through the CRT. These lines 
contain the date, run and point number, the PRT reading 
from the microvolt potentiometer with the last digit 
repeated, the reading of the pressure gage, the barometer 
reading, the wire resistance for both long and short hot 
wire, the sums of lead and ballast resistance for both long 
and short hot wire sections, the time increment at which the 
voltage readings were taken, and lastly the voltage applied 
to the bridge. The remainder of the file contains the set of 
250 voltage readings across the Wheatstone bridge. 

4.1. Cell Temperature 

The PRT is a standard capsule thermometer. It has been 
calibrated by the NBS temperature section. The voltage 
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FlCURE 4. Block diagram of temperature measurement and cryostat control circuits. 
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FIGURE 6. The minicomputer. 

read across the 10 $2 standard resistor is used to adjust the 
thermometer current so that a current of 1 ma is flowing in 
the PRT circuit. A second voltage reading across the PRT 
proper allows its resistance to be calculated. The second 
voltage reading is the only value that is keyed into the CRT. 
The error in the cell temperature as established from a 
calibration of the microvolt potentiometer used is 0.001 K 
at 80 K increasing to 0.007 K at 292 K. 

4.2. Cell Pressure 

The cell pressure is read with a commercial high preci- 
sion steel bourdon tube. The unit was selected because it is 
"compatible" with high pressure oxygen. The unit is cali- 
brated by the vendor with the calibration traceable to NBS. 
A digital readout is provided through an optical sensor, 
with the units of the readout in "counts." At the maximum 
pressure 96000 counts correspond to approximately 68 
MPa. At the higher pressures the bourdon tube displays 
hysteresis under loading as a function of time. It is this 
hysteresis that gives rise to the stated uncertainty of ± 0.03 
MPa. The pressure gage calibration was represented by a 
low order polynomial for use in the data reduction program. 
The error attributable to this curve fit is well within the 
stated error of the gage. Vapor pressure checks of the gage 
at approximately 1 MPa with both oxygen [26] and propane 
[27] confirm the calibration at the lower end of the pressure 
scale. 

4.3. Resistances 

Direct resistance measurements are made while balanc- 
ing the Wheatstone bridge. The resistances of the hot wires, 
of the leads, and of the ballast resistors are measured with 
the DVM set at the optimum gain. The resistance measure- 
ments are made using the four terminal method; the voltage 
across a standard 100 Q resistor defines the current in the 
circuit and the voltage reading across the unknown defines 
the resistance in question. In actual practice each voltage 



reading is the average of 100 individual readings. Each hot 
wire section measurement must be corrected for leads in- 
side the high pressure cell which are short pieces of copper 
and steel wire. The corrections were calculated from wire 
dimensions and tabulated resistivities [28], They are esti- 
mated to be 0.22 fl for the long hot wire section and 0.65 Q 
for the short hot wire section at room temperature. The cor- 
rections are handled in the data reduction program and are 
further adjusted for changes in bath temperature. 

In order to obtain the temperature increase of the plati- 
num wires from the corresponding resistance increase, we 
need to know the variation of resistance with temperature 
for both wires. It has been shown in the past [8, 12, 13, 29] 
that an in situ calibration of the wires is desirable and also 
that the resistances per unit length of both wires should be 
the same to within about 2 percent. 

The wire resistances measured at essentially zero applied 
power in the balancing of the bridge together with the cell 
temperatures as determined from the platinum resistance 
thermometer are taken as the in situ calibration of the 
wires. During measurements on oxygen [26] some 1800 
values were collected for each wire in the temperature range 
76-320 K with pressures from atmospheric to about 70 
MPa. The resistance relation for each wire is represented by 
an analytical function of the type 



R(t) = A + BT + Cr + DP 



(4) 



where T is the temperature in kelvin and P the pressure 
gage reading. The pressure dependence is small but statis- 
tically significant and reflects the fact that the calibration 
measurements are made with a small applied voltage of 0.1 
volts. Coefficients for eq (4) were determined in two ranges 
of temperature with some overlap in the range of the fits as 
shown in table 1. The standard deviation of the resistance 
measurements as well as the equivalent temperature errors 
are shown in table 1. The long wire has a length of 10.453 
cm at room temperature, the short wire one of 5.143 cm. 
Both wires have a nominal diameter of 0.001270 ± 
0.000001 cm, thus the radius a in eq (1) is 0.000625 cm. If 
we use these values and the measured resistances we can 
calculate a resistivity of 10.07 X 10" 6 fi-cm at 273.15 K com- 
pared to a best value of 9.60 X 10" 6 Q-cm for an annealed 
high purity specimen [28]. The difference is ascribed to the 
hard drawn condition of our wires. The a of these wires, 
defined as (i?(373.15)-i?(273.15)/(100. i?(273.15)) is 
0.0037944 which is lower than the value 0.003925 required 
for use in an annealed platinum resistance thermometer 
[26]. Finally, we evaluate 






and 



a -MEL 

os -m 



(5) 



463 



Table I. Calibration Constants of Wires. 
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29.3259 
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and compare a L and a s in the experimental range. Table 2 
shows the percent difference between o L and a s as a func- 
tion of T. It is apparent that for the lower temperatures the 
correction proposed by Kestin and Wakeham [8] would ap- 
ply except that in this apparatus independent calibrations 
for both long and short hot wire are used. 

4.4. Time Increment 

The timing of the experiment is based on the internal 
clock of the minicomputer, a temperature stabilized quartz 
crystal oscillator that divides a second into roughly 20,000 
parts. Available to the programmer is a real time clock 
which is incremented in milliseconds. During a run this 
clock is first zeroed, line 620 in appendix I, next 250 voltage 
readings are taken across the bridge, and finally the clock is 
read, line 670 in appendix I. The gain of the DVM is set 
before starting the voltage readings. There are 250 steps in 
the program loop, where each step includes reading the 
DVM, storing the result, and time for the computer to ex- 
ecute the step. The time required to execute the program 
loop is 755 ms, which means that the DVM is read every 
3.02 milliseconds. Separate timing runs establish the time 
necessary to read the clock as about 1 ms. This means the 
time increment could be as small as 3.012 ms. Actual calcu- 
lations using the data reduction program show that for this 



variation in time increment the thermal conductivity results 
are affected by less than 1 part in 20,000. 

The time increment between readings of the DVM can be 
varied by inserting a delay of several ms into the program 
loop. Increments of 5, 10, and 15 ms were tried. In general, 
the use of larger time increments simply reduces the num- 
ber of useful points for the thermal conductivity analysis 
since for nearly all fluid conditions the onset of convection 
occurs at around one second in this apparatus. 

4.5. Power Level 

The ideal analysis stipulates that the heat flux q, applied 
to the wire remains constant during a given run. Under ex- 
perimental conditions q is nearly ideal at room tempera- 
ture, i.e., it varies by about two parts in 10000. However, at 
the very lowest temperatures the q will vary from beginning 
to end of a run by values up to several percent. To see why 
this is so, and to find the correction that has to be applied 
consider the following. 

The power supply used provides a constant voltage at its 
terminals to a circuit which consists of a 100 12 standard 
resistor in series with the Wheatstone bridge. The effective 
bridge resistance is nominally 100 A but changes as the run 
progresses. The voltage of the power supply is preset, it is 
known to an accuracy of ± 200 /xV. The set value is entered 
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through the key board, just prior to switching the power to 
the bridge. After the run is completed, but before the power 
is switched back to the dummy resistor, the voltage across 
the standard resistor, channel in figure 2, is read. The dif- 
ference between this voltage and the preset one is the 
voltage applied to the bridge a short time after completion 
of the run. The difference is one of the items placed on 
magnetic tape, and it is used in the data reduction program 
to calculate the power applied to the bridge. The instan- 
taneous heat flux is determined from the equation. 



q (t) = pR/e = [V(t)i(R 3 (t) + R 4 (t))Y (Rrit) + 



(6) 



where V (t ) is the voltage applied to the bridge, and R 3 and 
R 4 are defined in figure 2. The voltage V(t), both hot wire 
legs of the bridge R 3 and i? 4 , and the hot wire resistances 
themselves change as a function of time. At room tempera- 
ture the drop in the current is almost exactly offset by the 
rise in the hot wire resistances. At low temperatures the 
fractional change in the hot wire resistances predominates. 
As an extreme example consider a run in liquid propane at 
say 110 K. For a 3 K rise in wire temperature the long hot 
wire rises by 1 12 from an initial value of 30 12, or 3.3 per- 
cent. This rise is only partially offset by the 0.7 percent drop 



in the hot wire current leaving a net rise in q of 2.6 percent 
from t = to t = 755 ms. The range of valid times for a run 
is generally t x = 150 ms, t 2 = 755 ms. Over this range of 
times the variation in q is less than the extreme case 
because a large part, 70 percent or so, of the resistance 
change of the hot wires occurs during the times not valid for 
the experiment. Nevertheless, the remaining variation in q 
is significant and we apply a correction to account for it. 
We select, a time in the middle of the valid time interval 
*<» = (*i + t 2 )l2 and correct the experimental AT's by the 
ratio of the instantaneous q to the q at the median time 



ATXfU, = AT(t) raw * q(t t )lq(t a ) . 



(7) 



V{t) varies by about one part in 1000 over a run, and some 
calculation verified that V(t) varies as the (n(t), and further 
that the value measured was actually achieved at the req- 
uisite time after the end of the run. 

4.6. Bridge Voltages 

The voltages measured across the bridge, actually a meas- 
ure of the bridge unbalance, are the basic data in this ex- 
periment. A typical set of 250 of these voltages for point No. 
9044 is shown in appendix II and is plotted in figure 7. The 
general shape of this curve is logarithmic, as expected. 




FIGURE 7. Typical bridge voltages versus time, point 9044, helium gas. 
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Noise levels in these readings are evident. Considerable ef- 
fort went into reducing the noise level, and also into tracing 
the source. The DVM, when shorted, exhibits an error band 
of ± 6 /tV. When a 100 fi resistor is placed across the DVM, 
no applied power yet, the band increases to ± 17 yN. With 
the experimental system connected to the DVM the error 
band becomes ± 28 fiY which is roughly equivalent to 
± .026 K in the AT. Only a small part, some 20 or 30 per- 
cent, is due to AC pickup. Use of integrating or averaging 
DVM's will reduce the noise level at the expense, however, 
of a smaller number of measurements. 



5. Data Reduction 

A listing of the data reduction program used for helium is 
given in appendix III. Experimental values are read from 
the magnetic tape. Initial calculations and corrections in- 
clude corrections for the bias of the DVM, corrections for 
the leads inside the high pressure cell, and an estimation of 
the instantaneous voltage applied to the bridge. The AT W 
values are calculated from the voltages applied to the 
bridge, the offset or measured voltages across the bridge 
and the two hot wire calibrations. The bridge equation is 
solved repeatedly until the increases in the hot wire 
resistances together with the voltage applied to the bridge 



yield an offset voltage equal to the measured one. The 
various corrections to the AT W are applied to obtain the AT, 
and then the AT's and £n(*)'s are passed to a least squares 
regression routine. Microfilm plots of both AT and the 
deviation of each measured point from the regression line 
versus tn(t) are set up. Typical examples again for point No. 
9044 are shown in figures 8a and 8b. In the final steps of the 
data reduction program an equation of state is used to find 
the applicable density, the experimental thermal conductiv- 
ity is compared to any previous correlation that might exist, 
and the results are printed. 

A certain amount of judgement enters into the data 
reduction process particularly in the selection of t t and t v 
the range of times over which the least squares straight line 
regression is to be evaluated. The criterion for a valid result 
is that the AT values as plotted in the example, figure 8a, 
shall form a straight line versus the tn(t). The deviation plot 
as shown in figure 8b is used to help determine the valid 
range of the least squares fit. For most experimental condi- 
tions t x = 150 ms and t 2 = 755 ms. In both figures, 8a and 
8b, the points plotted begin at index 1 1 or t = 33 ms. The 
first 30 or so points are seen to fall somewhat below the 
straight line. By inspection the fit for point 9044 could have 
included an index as low as 41 or t x - 124 ms. In figure 8b 
the standard deviation for the regression fit of point 9044 is 
roughly 0.2 percent, or at the 3a level 0.65 percent. 




4. a 5.6 
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FIGURE 8a. Typical temperature rises versus the logarithm of time, point 9044, helium gas. 
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FIGURE 8b. Typical deviations of experimental temperature rises from the calculated straight 
line versus the logarithm of time, point 9044, helium gas. 



By implication, it is equally important to recognize and 
interpret nonstandard results obtained, so that invalid 
results can be identified and rejected. A description of 
several of the more interesting malfunctions follows. In 
figure 9a we see that the AT values are nearly constant 
toward the end of the run. This state of affairs depending 
on the state of the fluid sample is interpreted as (a) the ap- 
paratus is operating in the steady state mode, or (b) as the 
onset of convection in the run. 

If a similar run at higher fluid densities is taken to much 
longer experimental times the plot is as shown in figure 9b. 
Here the effect of convection is clearly seen as a rapid cool- 
ing of the hot wires after the maximum AT is passed. In 
principle we could get a valid result for the transient hot 
wire by selecting t Y = 690 ms and t 2 = 1353 ms. However, 
the number of AT values useful for the straight line regres- 
sion has now dropped to 50. In practice it turns out that the 
applied power level, as seen in the resulting AT s of 10 to 12 
K, was simply too large to be useful. The plot of figure 9c 
shows the effect of gas circulation within the cell. The 
precise cause is not known. Convection combined with ther- 
mal oscillations when the capillary cools are suspected 
reasons. From the plot it is clear that the wires experience a 
non-uniform temperature field and are subject to abrupt 
segments of convection cooling. Measurements with a plot 
of this type are rejected. 



In figure 9d the DVM saturated toward the end of the 
run. By excluding the points for which the voltage readings 
are constant, and provided there are sufficient measure- 
ments to extract a straight line fit, a measurement of this 
type can be used. At the very lowest pressures (densities) in 
the vapor phase at temperatures much below the critical 
temperature some plots look like the one shown in figure 9e. 
The physical interpretation of this run is that we see a 
Knudsen effect at short times and an outer boundary effect 
at the longer times. No valid range for the hot wire instru- 
ment exists and these results are rejected. In the region 
near the critical point the onset of convection occurs much 
more rapidly. Measurements in this fluid region have to be 
made at lower power levels and at shorter times. For the 
measurements on oxygen experimental times of /, = 33 ms 
and t 2 = 500 ms were found to be appropriate. 

A second criterion which is used in accepting or rejecting 
experimental measurements is by considering the change in 
the measured thermal conductivity as the applied power, q, 
is varied. Examples are given in the sections on nitrogen 
and helium. In a general way this criterion suggests that 
measurements in which the final AT s are larger than about 
5 K are suspect. This includes most of the critical region 
where the results have to be evaluated quite carefully, 
almost on an individual point by point basis. 
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FlCURE 9a. Plots of experimental malfunctions. Point 29132, argon gas. 
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FlCURE 9b. Point 5027, nitrogen gas. 
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FIGURE 9c. Point 29010, argon gas. 




4.8 5.6 

ln( t ) . t m ms 



FIGURE 9d. Point 31022, propane liquid. 
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Figure 9e. Point 28003, oxygen vapor. 
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Figure 9f. Point 5035, nitrogen gas. 
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6. Performance Verification 

6.1. Nitrogen Results 

One of the gases selected for apparatus testing was 
nitrogen. Typical measurements on nitrogen are presented 
in table 3. There are two different sets of measurements. 
Run 4 is an abbreviated isotherm to check the functioning 
of the system as the pressure is varied. In contrast, all 
measurements of run 5 are taken at a single density while 
several of the pertinent apparatus parameters are varied. 

Run 4, the pseudo-isotherm, was measured at tempera- 
tures between 297 and 301 K with pressures up to 68 MPa. 
At each pressure (density) level a minimum of four different 
power settings were used. The results, a total of 50 points at 
varying temperatures are shown in table 3. For easy com- 
parison with the results of others and using a fixed dKldY of 
0.00063 W/m«K 2 the results were shifted at the experimen- 
tal densities to an even temperature of 300.65 K as shown in 
the next to last column of table 3. The results so adjusted 
are represented with a curve fit of the type used by Clifford, 
et al. [31] 



X = A + Bq + Cq 2 + Dq 3 



(8) 



The coefficients were determined as A = 0.02550, B — 
0.000112, C = 0.310567 X 10" 4 , andD = 0.261641 X 10" 5 
with X in W/m«K and q in mol/L. The standard deviation 
for the fit is 0.00022 W/m»K which translates to roughly 0.5 
percent in the middle of the density range. The experimen- 
tal measurements as adjusted to 300.65 K and the curve fit 
are shown in figure 10. Taking the curve fit of eq (8) as the 
base, the present results and those of others [31,32,33] are 
compared in figure 11. When all results are considered 
together they span a range of 3 percent. The differences 
between the present results and those of Clifford, et al. [31] 
who use a transient hot wire system range between 1.3 and 
2.3 percent, the differences to the more recent meas- 
urements of Assael, et al. [32] range from 1.1 to 1.7 percent. 
Differences between the present results and a published 
wide-range correlation for nitrogen [33] lie between 0.6 and 
4.5 percent. Reference [33] serves as the source of the equa- 
tion of state for nitrogen. 

Run 5 was measured at the same nominal pressure, about 
16.3 MPa, with experimental temperatures varying from 
297 to 309 K. For easy comparison these results are referred 
to a common temperature, 300.65 K, as indicated above. 
However, they are further adjusted to a common density of 
6.3 mol/L using eq (8). The values so adjusted are shown in 
the last column of table 3, and are plotted in figure 12 ver- 
sus the power applied, q. The first 15 points, shown as 
circles in figure 12, are repetitive measurements to establish 
the base of comparison. The mean of these 15 points is the 



straight line in the center of figure 12 while the lines above 
and below represent changes of one percent. 

We indicated previously that the time increment can be 
varied. A series of such measurements was made on 
nitrogen gas with time delays of 5 and 10 ms. The runs are 
shown in table 3 as points 5016 — 5028. The actual time in- 
crement between voltage measurements is 9.02 ms for runs 
5016—5020 and 14.02 ms for runs 5021—5028 compared to 
a normal time increment of 3.02 ms. A plot of point 5027 is 
given in figure 9b, where the straight line segment indicates 
the range of times, t x = 154.22 ms and t 2 = 757.08 ms, used 
to evaluate the result given in table 3. If proper times t x and 
t 2 are selected, then runs of this type can yield valid results. 
However, the associated regression statistics deteriorate 
because a smaller number of points is available to the least 
squares analysis. 

A second test involves larger than normal powers. The 
idea is to compare measurements made at one reference 
temperature with a very large power to measurements made 
at a higher reference temperature with a much smaller ap- 
plied power, an overlapping of isotherms so to speak. For 
example a run at T re f=295 K and a power level to produce 
A7"s of around 15 K or a T exp of 310 K could be compared 
to a second run of T rtf = 307 K and a power level to pro- 
duce Ars ~ 4 K with a T T „ also of 310 K. The 
measurements taken are shown in table 3 as points 5031 — 
5035. The test fails because we can see a slight dependence 
on power level for the high powers in figure 12. The reason 
the test fails is because the plot of AT versus (n(t) instead of 
being a straight line is curved at the very large ATs in- 
volved, as shown in figure 9f for point 5035. The test might 
become valid in the future, provided a reanalysis of the cor- 
rections to the experimental AT W is made which allows for 
sufficiently large AT. 

A final experimental test involves the gradient along the 
high pressure cell. The cell is fitted with a thermocouple 
that measures the temperature difference from the top to 
the bottom of the cell. Nearly all runs were measured with 
the top of the cell slightly warmer than the bottom, about 10 
to 15 mk. Test runs with the cell top hotter by 47 mk are 
shown as points 5036 — 5040 while similar measurements 
with the cell top colder by about 35 mk are given by points 
5041 — 5045. The results show that the thermal conductiv- 
ities measured under these conditions vary by roughly 0.2 
percent, which in view of the other uncertainties is 
negligible. 

6.2. Helium Results 

The rare gas selected initially was helium. Results ob- 
tained subsequently on argon have already been published 
[30]. For helium a single pseudo isotherm was measured at a 
nominal temperature of 307 K with pressures up to 68 MPa. 
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Table 3. Thermal Conductivity of Nitrogen Near 300 K. 



Run Pt. 



4001 
4002 
4003 
4004 
4005 
4006 
4007 
4008 
4009 
4010 
4011 
4012 
4013 
4014 
4015 
4016 
4017 
4018 
4019 
4020 
4021 
4022 
4023 
4024 
4025 
4026 
4027 
4028 
4029 
4030 
4031 
4032 
4033 
4034 
4035 
4036 
4037 
4038 
4039 
4040 
4041 
4042 
4043 
4044 
4045 
4046 
4047 
4048 
4049 
4050 

5001 
5002 
5003 
5004 
5005 
5006 
5007 
5008 
5009 
5010 
5011 
5012 
5013 
5014 
5015 



Pressure 
MPa 

69.123 

69.124 

69.124 

69.123 

69.125 

69.125 

61.132 

61.133 

61.128 

61.131 

54.586 

54.588 

54.588 

54.588 

47.903 

47.903 

47.905 

47.905 

40.607 

40.608 

40.609 

40.609 

33.596 

33.596 

33.596 

33.596 

26.643 

26.642 

26.643 

26.642 

19.700 

19.700 

19.700 

19.700 

19.700 

12.708 

12.708 

12.708 

12.708 

5.681 

5.681 

5.681 

5.681 

5.681 

5.681 

1.430 

1.430 

1.429 

1.429 

1.429 

16.327 
16.327 
16.327 
16.327 
16.327 
16.328 
16.328 
16.329 
16.329 
16.329 
16.329 
16.329 
16.329 
16.329 
16.329 



Temperature 
K 

297 .004 
297.439 
298.096 
298.759 
299.583 
300.322 
298.146 
298.905 
299.709 
300.741 
298.026 
298.703 
299.475 
300.334 
300.778 
299.771 
298.863 
298.042 
297.907 
298.705 
299.584 
300.704 
298.229 
299.023 
300.028 
301.176 
297.741 
298.530 
299.522 
300.611 
300.643 
299.992 
298.934 
297.950 
297.532 
297.418 
298.319 
299.457 
300.749 
301.332 
300.556 
299.102 
297.843 
299.714 
298.464 
299.675 
298.299 
297.129 
297.703 
298.984 

297.649 
298.523 
299.492 
300.657 
301.440 
301.373 
300.632 
299.496 
298.426 
297.614 
297.512 
298.293 
299.270 
300.472 
301.184 



Density 
mol/L 

17.5101 

17.4916 

17.4635 

17.4351 

17.4003 

17.3691 

16.4493 

16.4166 

16.3813 

16.3376 

15.5156 

15.4863 

15.4528 

15.4156 

14.3096 

14.3531 

14.3927 

14.4285 

13.0618 

13.0277 

12.9903 

12.9427 

11.4980 

11.4655 

11.4247 

11.3782 

9.7001 

9.6705 

9.6340 

9.5939 

7.4708 

7.4902 

7.5221 

7.5521 

7.5651 

5.0731 

5.0546 

5.03i4 

5.0054 

2.2732 

2.2798 

2.2923 



2.3032 
2.2870 
2.2975 
.5751 
.5779 
.5799 
.5788 
.5761 

6.3985 
6.3757 
6.3507 
6.3209 
6.3010 
6.3030 
6.3218 
6.3511 
6.3787 
6.3999 
6.4028 
6.3824 
6.3571 
6.3263 
6.3083 



Power 
W/m 

.14023 
.20067 
.27180 
.35381 
.44662 
.55006 
.27174 
.35371 
.44650 
.54992 
.23468 
.31117 
.39862 
.49695 
.49683 
.39857 
.31119 
.23469 
.20044 
.27163 
.35361 
.44624 
.20061 
.27208 
.35423 
.44717 
.14033 
.20078 
.27210 
.35421 
.31180 
.27212 
.20073 
.14029 
.11414 
.09066 
.14029 
.20079 
.27209 
.23508 
.20076 
.14028 
.09064 
.16920 
.11412 
.14031 
.09064 
.05179 
.06983 
.11411 

.09069 
.14036 
.20084 
.27209 
.31171 
.31172 
.27203 
.20069 
.14025 
.09061 
.09071 
.14037 
.20090 
.27222 
.31195 



Thermal 

Conductivity 

W/m.K 

.06868 
.06880 
.06821 
.06829 
.06845 
.06831 
.06344 
.06379 
.06333 
.06347 
.05971 
.05962 
.05980 
.05977 
.05555 
.05551 
.05570 
.05555 
.05049 
.05078 
.05092 
.05104 
.04616 
.04600 
.04636 
.04633 
.04128 
.04135 
.04199 
.04154 
.03707 
.03693 
.03695 
.03678 
.03676 
.03223 
.03212 
.03216 
.03238 
.02837 
.02830 
.02816 
.02813 
.02821 
.02807 
.02605 
.02582 
.02618 
.02581 
.02608 

.03439 
.03453 
.03484 
.03489 
.03508 
.03485 
.03495 
.03467 
.03471 
.03459 
.03447 
.03453 
.03463 
.03468 
.03482 



STAT 



.034 
.021 
.013 
.009 
.007 
.004 
.012 
.009 
.006 
.005 
.014 
.009 
.007 
.006 
.004 
.006 
.009 
.014 
.016 
.010 
.007 
.005 
.015 
.010 
.007 
.005 
.023 
.013 
.009 
.006 
.006 
.008 
.012 
.019 
.026 
.035 
.017 
.011 
.008 
.010 
.013 
.021 
.042 
.016 
.030 
.024 
.041 
.097 
.064 
.030 

.032 
.019 
.011 
.007 
.006 
.006 
.006 
.011 
.018 
.035 
.032 
.018 
.011 
.007 
.006 



Adjusted Thermal Conductivity 

with T = 300.65 K 
with T = 300.65 K and p = 6.3 mol/L 
W/m.K W/m»K 

.06891 
.06900 
.06837 
.06841 
.06852 
.06833 
.06360 
.06390 
.06339 
.06346 
.05988 
.05974 
.05987 
.05979 
.05554 
.05557 
.05581 
.05571 
.05066 
.05090 
.05099 
.05104 
.04631 
.04610 
.04640 
.04630 
.04146 
.04148 
.04206 
.04154 
.03707 
.03697 
.03706 
.03695 
.03696 
.03243 
.03227 
.03224 
.03237 
.02833 
.02831 
.02826 
.02831 
.02827 
.02821 
.02611 
.02597 
.02640 
.02600 
.02618 



.03458 
.03466 
.03491 
.03489 
.03503 
.03480 
.03495 
.03474 
.03485 
.03478 
.03467 
.03468 
.03472 
.03469 
.03479 



.03440 
.03452 
.03482 
.03485 
.03503 
.03479 
.03491 
.03465 
.03471 
.03460 
.03448 
.03453 
.03462 
.03464 
.03477 
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Table 3 


. (Continued) 
Thermal 




tun Pt. 


Pressure 


Temperature 


Density 


Power 


Conductivity 


STAT 




MPa 


K 


raol/L 


W/m 


W/m.K 




5016 


16.330 


297.460 


6.4043 


.09069 


.03505 


.060 


5017 


16.330 


298.264 


6.3834 


.14035 


.03459 


.033 


5018 


16.330 


299.300 


6.3566 


.20085 


.03481 


.019 


5019 


16.330 


300.512 


6.3255 


.27213 


.03476 


.012 


5020 


16.330 


301.228 


6.3074 


.31179 


.03470 


.009 


5021 


16.329 


296.862 


6.4199 


.05183 


.03397 


.186 


5022 


16.329 


297.554 


6.4017 


.09073 


.03459 


.081 


5023 


16.329 


298.358 


6.3807 


.14043 


.03436 


.041 


5024 


16.329 


299.454 


6.3524 


.20091 


.03447 


.024 


5025 


16.329 


300.573 


6.3235 


.27225 


.03468 


.014 


5026 


16.329 


302.013 


6.2871 


.35420 


.03497 


.011 


5027 


16.329 


303.570 


6.2483 


.44706 


.03505 


.009 


5028 


16.329 


305.150 


6.2095 


.55073 


.03486 


.010 


5031 


16.328 


301.848 


6.2910 


.35473 


.03465 


.008 


5032 


16.328 


303.341 


6.2537 


.44782 


.03481 


.006 


5033 


16.328 


305.077 


6.2109 


.55174 


.03514 


.005 


5034 


16.328 


306.934 


6.1660 


.66648 


.03509 


.004 


5035 


16.327 


308.967 


6.1175 


.79195 


.03528 


.003 


5036 


16.333 


297.573 


6.4025 


.09074 


.03478 


.037 


5037 


16.333 


298.377 


6.3815 


.14047 


.03471 


.019 


5038 


16.333 


299.450 


6.3537 


.20096 


.03479 


.011 


5039 


16.333 


300.685 


6.3222 


.27231 


.03488 


.007 


5040 


16.333 


301.270 


6.3073 


.31209 


.03483 


.006 


5041 


16.337 


297.448 


6.4073 


.09066 


.03435 


.034 


5042 


16.337 


298.207 


6.3874 


.14029 


.03464 


.019 


5043 


16.337 


299.265 


6.3600 


.20074 


.03469 


.011 


5044 


16.338 


300.475 


6.3293 


.27202 


.03477 


.007 


5045 


16.338 


301.173 


6.3115 


.31171 


.03483 


.006 



Adjusted Thermal Conductivity 

with T - 300.65 K 
with T - 300.65 K and P - 6.3 mol/L 
W/m.K W/m-K 



.03525 
.03474 
.03490 
.03477 
.03466 
.03421 
.03479 
.03450 
.03455 
.03468 
.03488 
.03487 
.03458 
.03457 
.03464 
.03486 
.03469 
.03476 
.03497 
.03485 
.03487 
.03488 
.03479 
.03455 
.03479 
.03478 
.03478 
.03480 



.03506 
.03459 
.03480 
.03472 
.03465 
.03399 
.03460 
.03435 
.03445 
.03464 
.03490 
.03496 
.03474 
.03459 
.03472 
.03502 
.03493 
.03509 
.03478 
.03470 
.03477 
.03484 
.03478 
.03435 
.03463 
.03467 
.03473 
.03478 
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FIGURE 10. Nitrogen measurements shifted to 300.65 K with curve fit, eq (8). 
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FIGURE 11. Comparison of the present results with those of others for nitrogen at 300.65 K. 

Base line eq (8), A this paper, □ Clifford, et al. [31], O Assael, et al. [32], Hanley, et al. 

[33] correlation. 



At each pressure (density) level two different power settings 
were used, while at three selected densities a larger range of 
power was covered. The results, a total of 77 points at vary- 
ing temperatures are shown in table 4. For easy comparison 
with the results of others and using a fixed dX/dT of 
0.00035 W/m»K 2 the results were shifted at the experimen- 
tal densities to a temperature of 300.65 K as shown in the 
last column of table 4. The results so adjusted are repre- 
sented with a curve fit of the same type as used on argon 
earlier. 



X = A + Bq + C {e 013 > - 1.0} 



(9) 



The coefficients were determined as A = 0.153099, B = 
0.108257 X 10" 2 , and C = 0.613662 X 10" 3 with X in 
W/m»K and q in mol/L. The standard deviation for the fit is 
0.00049 W/m»K which translates to roughly 0.6 percent at 
the 95 percent confidence level or 0.9 percent at the 99 per- 
cent confidence level. In other words, of the 77 measure- 
ments 71 lie within ± 0.5 percent of the fitted curve. The 
experimental measurements as adjusted to 300.65 K and the 
curve fit are shown in figure 13. The densities at which a 
larger variety of power settings was measured are easily 



noted in figure 13 as 5.7, 11.5, and 19.6 mol/L. Taking the 
curve fit of eq (9) as the base, the present results and those 
of others [7,34,35,36,37] are compared in figure 14. These 
results span a range of ±2 percent. The deviations between 
the present results and those of Le Neindre, et al. [34], 
Johannin, et al. [35], and Ho and Leidenfrost [36] are a 
nominal 1 percent or less. The differences between the pres 
ent results and those of Kestin, et al. [7], who also use a tran 
sient hot wire system, range between 1.3 and 1.9 percent 
while the differences between the present results and a pub 
lished wide-range correlation for helium [37] also lie be 
tween 1 and 2 percent. Differences between the present 
results and the most recent ones of Assael, et al. [10], who 
also use a transient hot wire system, range from 1.0 to 1.7 
percent. The source used as equation of state for obtaining 
the densities of helium is [37]. 

Finally, we apply a test that has not been applied exten- 
sively to other transient hot wire systems. In analogy to a 
flat plate steady state thermal conductivity apparatus [38] 
we insist that the measured thermal conductivity in the 
present apparatus should be independent of the applied 
power. Figure 15 shows the measurements for helium gas at 
19.6, 1 1.5, and 5.7 mol/L plotted as a function of the applied 
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FIGURE 12. Measured thermal conductivity as a function of the applied power. Apparatus 

tests with nitrogen gas. 

O normal measurements, A cell top cold, V cell top hot 

□ 5 ms delay, time increment 9.02 ms 

10 ms delay, time increment 14.02 ms 

* larger than normal powers. 



power. Whereas q changes by a factor of nearly five no 
dependence of the measured conductivity on power level is 
discernible. For reference, the values of the curve fit (eq (9)) 
plus or minus 0.5 percent are shown as three straight lines 
for each density in figure 15. The error bars for each point 
are those originating from the regression straight line fit for 
each point, see column STAT in table 4. 

6.3. Cell Alinement 



applied power, q. If the cell is out of alinement the 
measured thermal conductivity increases. The increase is 
more pronounced at the higher power levels. The proper 
setting for the alinement screws occurs when the measured 
results are most nearly independent of q as illustrated in 
figure 15. The final settings chosen were nearly identical to 
those first achieved during construction of the apparatus. In 
fact, the measurements on helium, Run No. 9, were taken 
four months before the series on alinement, Run No. 19. 



During the construction of the apparatus we sought to 
keep the cell in vertical alinement to insure that convection 
in the cell would be delayed as long as possible. Later we 
sought to verify the alinement by a series of experimental 
measurements, some 162 points with oxygen gas as a sample 
at a pressure near 10 MPa and a temperature of around 240 
K. The position of each cell alinement screw was varied in 
turn by a measured amount. For each setting of the screws a 
series of seven to nine measurements were taken varying the 



6.4. Other Tests 

Tests were run on nitrogen gas at room temperature with 
only a single hot wire in the bridge. The circuit was rewired 
so that either R A or R 3t see figure 2, were replaced by yet 
another 100 standard resistor. The data reduction pro- 
gram had to be modified accordingly. The results were eval- 
uated only on the minicomputer, and are relative to similar 
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Table 4. Thermal Conductivity of Helium Near 300 K. 



Run Pt. 



9001 
9002 
9003 
9004 
9005 
9006 
9007 
9008 
9009 
9010 
9011 
9012 
9013 
9014 
9015 
9016 
9017 
9018 
9019 
9020 
9021 
9022 
9023 
9024 
9025 
9026 
9027 
9028 
9029 
9030 
9031 
9032 
9033 
9034 
9035 
9036 
9037 
9038 
9039 
9040 
9041 
9042 
9043 
9044 
9045 
9045 
9046 
9047 
9048 
9049 
9050 
9051 
9052 
9053 
9054 
9055 
9056 
9057 
9058 
9059 
9060 
9061 
9062 
9063 
9064 
9065 



Pressure 
MPa 

63.980 
63.978 
63.978 
63.979 
63.979 
63.979 
61.905 
61.905 
59.692 
59.693 
57.621 
57.622 
55.588 
55.588 
53.338 
53.338 
51.292 
51.293 
49 .069 
49.070 
46.963 
46.964 
44.676 
44.678 
42.449 
42.449 
40.222 
40.222 
38.053 
38.053 
35.850 
35.851 
33.592 
33.592 
33.592 
33.592 
33.593 
33.594 
33.594 
33.594 
33.594 
33.595 
33.595 
33.595 
33.595 
31.497 
31.497 
28.968 
28.968 
26.634 
26.634 
24.532 
24.532 
22.396 
22.395 
20.255 
20.255 
17.991 
17.913 
15.627 
15.627 
15.626 
15.626 
15.626 
13.436 
13.436 



Temperature 
K 

306.154 
306.459 
306.788 
307.287 
307.741 
308.241 
306.849 
308.397 
306.827 
307.794 
306.382 
307.822 
306.292 
307.767 
306.324 
307.769 
306.429 
307.763 
306.377 
307.784 
306.311 
307.749 
306.361 
307.731 
306.396 
307.874 
306.328 
307.816 
306.321 
307.861 
306.292 
307 .886 
305.206 
305.508 
305.853 
306.306 
306.741 
307.256 
307.857 
305.799 
306.338 
306.824 
307.313 
307.790 
308.438 
306.173 
307.807 
306.381 
307.941 
306.344 
307.965 
306.354 
307.963 
306.379 
308.019 
306.511 
308.186 
306.347 
308.066 
305.967 
306.425 
306.897 
307.578 
308.170 
305.988 
307.519 



Density 
mol/L 



.6780 
.6629 
.6469 
.6229 
.6010 
.5769 
.1335 
.0601 
.5826 
.5380 
8.0798 
8.0144 
5633 
7.4977 
6.9781 
6.9154 
6.4358 
6.3793 
5.8453 
5.7876 
5.2786 
5.2211 
4.6495 
4.5966 
4.0279 
3.9726 
3.4010 
3.3472 
2.7787 
2.7251 
2.1378 
2.0847 
1.5040 
1.4943 
1.4833 
1.4688 
1.4552 
1.4389 
1.4198 
1.4856 
1.4684 
1.4531 
1.4375 
1.4223 
1.4018 
0.8425 
0.7931 
0.0629 
0.0184 
9.3369 
9.2936 
8.6709 
8.6306 
7.9829 
7.9445 
7.2793 
7.2434 
6.5288 
6.4692 
5.7333 
5.7254 
5.7170 
5.7052 
5.6951 
4.9748 
4.9516 



Power 
W/m 

.36457 
.46029 
.56726 
.68489 
.81392 
.95386 
.56688 
.95337 
.56708 
.81386 
.46038 
.81384 
.46046 
.81391 
.46030 
.81397 
.46023 
.81388 
.46028 
.81385 
.46031 
.81389 
.46023 
.81401 
.46025 
.81371 
.46013 
.81378 
.46027 
.81376 
.46029 
.81381 
.20675 
.28015 
.36455 
.46028 
.56704 
.68485 
.81389 
.36407 
.45965 
.56620 
.68391 
.81285 
.95255 
.46033 
.81402 
.46036 
.81396 
.46035 
.81397 
.46038 
.81417 
.46046 
.81418 
.46102 
.81536 
.46009 
.81378 
.36458 
.46017 
.56700 
.68465 
.81363 
•36453 
.68483 



Thermal 

Conductivity 

W/m.K 

.18386 
.18445 
.18336 
.18479 
.18424 
.18427 
.18268 
.18285 
.18132 
.18250 
.17957 
.18141 
.17811 
.17947 
.17822 
.17892 
.17730 
.17749 
.17662 
.17682 
•17652 
.17577 
.17437 
.17554 
.17243 
.17387 
.17191 
•17336 
.17121 
.17216 
.17026 
.17134 
.16850 
.16954 
.16908 
.16863 
.17039 
.16933 
.17057 
.16952 
•16955 
.16952 
.16982 
.17030 
.17128 
.16840 
•16980 
.16811 
.16769 
•16695 
.16735 
.16588 
•16621 
•16482 
.16552 
•16399 
.16490 
.16298 
.16397 
•16262 
.16148 
.16220 
.16218 
.16305 
.16032 
.16199 



STAT 



.022 
.014 
.011 
.009 
.006 
.005 
.010 
.005 
.010 
.007 
.014 
.006 
.015 
.006 
.015 
.007 
.015 
.006 
.014 
.006 
.015 
.006 
.013 
.006 
.013 
.006 
.014 
.006 
.013 
.006 
.014 
.006 
.042 
.029 
.020 
.014 
.010 
.008 
.006 
.020 
.013 
.011 
.008 
.006 
.005 
.013 
.006 
.013 
.006 
.014 
.005 
.014 
.006 
.013 
.006 
.013 
.006 
.014 
.006 
.018 
.012 
.010 
.008 
.006 
.019 
.007 



Adjusted Thermal 

Conductivity 

with T = 300.65 K 

W/m.K 

.18193 
.18242 
.18121 
.18247 
.18176 
.18161 
.18051 
.18014 
.17916 
.18000 
.17756 
.17890 
.17614 
.17698 
.17623 
.17643 
.17528 
.17500 
.17462 
.17432 
.17454 
.17329 
.17237 
.17306 
.17042 
.17134 
.16992 
.17085 
.16923 
.16964 
.16829 
.16881 
.16691 
.16784 
.16726 
.16665 
.16826 
.16702 
.16805 
.16772 
.16756 
.16736 
.16749 
.16780 
.16855 
.16647 
.16730 
.16610 
.16514 
.16496 
.16479 
.16388 
.16365 
.16281 
.16294 
.16194 
.16226 
.16099 
.16137 
.16076 
.15946 
.16001 
.15976 
.16042 
.15841 
.15959 
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Table 4. (Continued) 



Run Pt. 



9066 
9067 
9068 
9069 
9070 
9071 
9072 
9073 
9074 
9075 
9076 















Adjusted Thermal 










Thermal 




Conductivity 


Pressure 


Temperature 


Density 


Power 


Conductivity 


STAT 


with 


T - 300.65 K 


MPa 


K 


mol/L 


W/m 


W/m«K 






W/m#K 


11.253 


305.990 


4.2056 


.36453 


.16059 


.019 




.15872 


11.253 


307.633 


4.1844 


.68454 


.16030 


.008 




.15786 


9.013 


305.969 


3.4013 


.36451 


.15889 


.019 




.15703 


9.013 


307.727 


3.3827 


.68445 


.15957 


.007 




.15709 


6.648 


305.860 


2.5362 


.36460 


.15773 


.012 




.15591 


6.648 


307.452 


2.5232 


.68476 


.15851 


.006 




.15613 


4.466 


305.982 


1.7199 


.36440 


.15684 


.016 




.15497 


4.466 


307.438 


1.7119 


.68485 


.15705 


.008 




.15467 


2.419 


306.028 


.9401 


.36441 


.15672 


.023 




.15484 


2.418 


307.670 


.9349 


.68458 


.15595 


.012 




.15349 


.415 


306.143 


.1627 


.36423 


.15477 


.222 




.15285 




Dens i t y ,mo I / I 
FIGURE 13. Helium measurements shifted to 300.65 K with curve fit, eq (9). 
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curvefit value 



r i.C4 



19.6 mol/L 



curvefit value 



curvefit value 



-o 



11.5 mol/L 



$==£ 



5.7 mol/L 



0.0 .2 .a .6 .8 1.0 

Heat Flux, W/m 
FlCURE 15. Measured thermal conductivity as a function of the applied power. Apparatus test 
with helium gas. 
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measurements made the same day with the full bridge. The 
apparent conductivities for the single wires increase above 
those of the full bridge as follows 



full bridge 1.000 

long hot wire only 1.048 

short hot wire only 1.083. 

Measurements of this type show the relative effect of the 
end corrections which would be required for the single wire, 
but which are presumed to be off setting in the bridge 
circuit. 

7. Errors 

Precision, reproducibility, internal consistency, and accu- 
racy are discussed here. 



pane seven different liquid isotherms can be represented 
with an overall deviation of 1.2 percent (2a). For oxygen, 13 
isotherms from the surface which includes gas, liquid, sub- 
critical vapor and the near critical region for an average 
deviation of 1.6 percent. 

Accuracy. The accuracy of the measurement can in prin- 
ciple be established from the measurements and certain 
theoretical considerations, i.e., for the rare gases the 
Eucken factor [7]. For argon [30] the 0.3 percent obtained in 
the low density extrapolation was considered a fortunate 
coincidence and an overall accuracy of ±1 percent was 
claimed. For the results on helium in table 4 the low density 
extrapolation yields a result which is 1.8 percent low. Accu- 
racy can also be estimated by comparison to the results of 
others. For the present apparatus these intercomparisons 
cover a wide range of fluids, and a wide range of pressures, 
densities and temperatures. Summarized these comparisons 



Precision. The measurement variation in the voltages 
across the Wheatstone bridge is ±28 /*V which corresponds 
to about 25 mK. The effective errors in the experimental 
AT values are 6 percent for the very small power levels, i.e., 
a final AT of about 0.5 K, to 0.4 percent for a final AT of 
about 6 K. An example of the distribution of the exper- 
imental ATs about the calculated straight line was shown 
in figure 8b. The regression program accounts for the 
spread in measurements by assigning a variance to the 
coefficients of the calculated straight line. The value STAT 
as printed in tables 3 and 4 is the fractional error in the 
slope of the straight line obtained from the variance in the 
slope (2a) divided by by the slope itself. The value STAT is 
the primary measure of the experimental precision because 
other contributions such as an error in q are quite negligi- 
ble. For final AT values of 3.5 K or more which includes 
about 70 percent of all measurements made to date the 
nominal precision (2a) is 0.6 percent or slightly less. 

Reproducibility was tested by measuring (a) repeatedly 
the same point at the same fluid conditions on a given day, 
(b) by measuring similar points at different power levels, 
and (c) by repeating different isotherms at different times 
some months apart. Examples of (a) are found in table 3, 
examples of (b) are given in the sections on nitrogen and 
helium. For (c), isotherms at 160 K, 200 K, and 240 K were 
repeated during the measurements on oxygen [26]. 

Internal consistency. When measurements are made over 
a large portion of the thermal conductivity surface, and if 
an analytical function is used to represent the surface, then 
it becomes possible to check the internal consistency of the 
set of measurements by intercomparing the deviations and 
deviation patterns for the various isotherms. Such checks 
are possible for both propane [27] and oxygen [26]. In pro- 



Nitrogen: this apparatus is 1.3 to 2.3 percent lower than 

other authors [31,32,33]. 
Helium: this apparatus is 1. to 1.8 percent lower than 

other authors [7,34,35,36,37]. 
Argon: this apparatus is 1. percent lower than other 

authors [30]. 
Oxygen: this apparatus is 1-2 percent lower than other 

authors [26]. 
Propane: this apparatus is higher at some densities, lower 

at other densities when compared to other 

authors [27]. 
Hydrogen: this apparatus is 1.5 percent lower than other 

authors [39]. 

8. Summary 

A new transient hot wire thermal conductivity apparatus 
for fluids has been developed. It differs from some of the 
other transient hot wire systems in that it uses compen- 
sating wires in a Wheatstone bridge arrangement. It differs 
from those transient hot wire systems that use a bridge in 
that the voltages developed across the bridge are measured 
directly. The apparatus is a very rapid measuring device. 
Sixty points, a full pseudo-isotherm, are possible in a day. 
The system is capable of measuring a wide range of fluid 
temperatures, densities, and pressures. Fluid states, such as 
dilute gas, dense gas, the near critical states, vapor at tem- 
peratures below critical, compressed liquid states, and 
metastable liquid states at densities below saturation, are 
accessible for many liquefied gases. A wide variety of simple 
fluids has been measured in the system so far. It appears 
that the precision of the new system is on the order of 0.6 
percent while its accuracy is close to 1.5 percent. 
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Appendix I. The Minicomputer Program 



10 REM PROGRAM TEST? RESULTS ONTO TAPE? REVISED 20 DEC 79 

20 DIM AL'2553jBI:2553 

30 LET 12=250 

40 REM WHEATSTONE BALANCE 

50 REM HOTWIRE SIDE ONLY 

60 PRINT " DIAL IN VOLTS SHOULD BE 0*1" 

70 PRINT "TOTAL R3"v "TOTAL R4"*'R3 POT SET'»'R4 POT SET" 

75 CALL 9>3* 

80 CALL 9,1y6 

90 FOR J= TO 7 

1.00 LET vl:j:]= 

105 CALL 5>l»128rJ>J*4 

110 FOR 1=1 TO 250 

120 CALL 6tAt.Il 

i so let v c j ::i =v i: j 3 + a i: :i: :i 
160 next i 
170 let ui:j:j:=vi::j:i/25o 

180 next j 

190 CALL 9»:l.y 

210 LET C=VC 03/100 

220 FOR 1=1 TO 7 

230 LET RCI3=VCI3/C 

240 NEXT I 

250 LET R3=RC13+RII23+RC33 

260 LET R4=RC43+RE:53+RC73 

300 print R3yR4?Ri::2::i.»Ri:5:] 

310 PRINT "GO AGAIN ? " 

320 INPUT KM. 

330 IF Kl= GOTO 80 

340 REM FULL BRIDGE 

350 LET VC63 = 

360 CALL 9*3»6 

370 CALL 9?:l. ,6 

380 CALL 5*lrl024>6>6»4 

390 FOR 1=1 TO 250 

400 CALL 6 r AIM 3 

410 let ui::6:i=vi:6:ifAci.i 

420 NEXT I 

430 let vi:6::i=mi:6:i/250 

440 CALL 9?1 j- 

450 CALL 9*3, 

*&0 PRINT "WHEATSTONE BALANCE =",VC63»" GO AGAIN ? ■ 

470 INPUT Kl 

480 IF Kl= GOTO 350 

1 9 R E M IN I T I A I... C N D I T 1 N S 

500 LET Tl= 

510 LET T2=296*06f (Ri::3:i~90* 58)/ ♦ 33 

520 LET AC 03-VE63 

530 PRINT 

540 PRINT " TIME-Z " v " TEMP-Z " r " R3-Z " r " R4~Z " r H V OFFS-Z " 

550 PRINT Tl,T2»RC33fftC43rAC 03 

48] 



560 PRINT 

570 REM DATA AQUISITION SECTION 

580 PRINT " SET HP POWER" 

590 INPUT VI 

600 CALL 9,3,6 

610 CALL 5,1,512,6,6,9 

620 CALL 1? 0, 

630 CALL 9,1,6 

640 FOR 1=1 TO 12 

650 CALL 6yAIIi::i 

660 NEXT I 

670 CALL 2,A2,A3,A4,A5 

680 FOR 1=1 TO 100 

690 CALL 7, 0»BIII.l 

700 NEXT I 

710 LET V2= 

720 FOR 1=1 TO 100 

730 LET y2=V2+BCI3 

740 NEXT I 

750 LET V2=V1~V2/10<> 

760 CALL 9,1, 

770 CALL 9,3, 

780 LET A4=A4+A5/1000 

790 LET A5=A4/.T2 

795 GOTO 900 

800 PRINT 

810 PRINT 'TOTAL RUN TIME", "TIME [NCR," 

820 PRINT A4,A5 

821 PRINT 

823 PRINT "DATE RUN/PT PRT - L PRESS BAR '?" 

824 INPUT 01,02,03,04,05,06 

825 LET S4=R3-RC33 

826 LET S5=R4-RC43 

827 PRINT 

828 CALL 22 

829 PRINT 01 5", "02? ","03? ","04 J", "05?", "06?"," 

830 PRINT RC33;","S4;","RC43J",»S5?","A5?","^;"," 
840 FOR 1=1 TO 49 

850 LET N=<I-1)*5 

870 NEXT 1 ! 1 ACN+1:] ? " ' " ACN+2:1 ? " » " ACN+3:i ' " ' " ACN+4J 5 " , " ACN +53 ? 

872 lli^^ 2 * 62 *'*'*^ 
874 END 
900 PRINT 

920 PRINT AC2463? "'" AC247:,? "' HAC248 ^"'"An249.1 J ","AC250_1 

930 PRINT "RECORD ON TAPE'?" 

940 INPUT Kl 

945 PRINT 

950 IF Kl= GOTO 800 

960 GOTO 874 
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Appendix II. Data for Point 9044, Helium Gas 



800116, 9044, 57515.7, .7, 4710.1, 
93.3094, 6.65579, 46.7144, 53.4657 
7.06954E-4, 5.15612E-4, 1.07951E-3 
2.00463E-3, 2.1614E-3, 
2.55064E-3, 2.65734E-3 
2.91842E-3, 2.96074E-3 
3.16639E-3, 3.16878E-3 
3.35952E-3, 3.38992E-3 
3.55206E-3, 3.56219E-3 
3.6671E-3, 3.6975E-3, : 
3.77559E-3, 3.84712E-3 
3.90017E-3, 3.90673E-3 
3.96335E-3, 3.99375E-3 
4.02296E-3, 4.09568E-3 
4.13562E-3, 4.14874E-3 
4.20119E-3, 4.19881E-3 
4.27332E-3, 4.27391E-3 
4.35677E-3, 4.34485E-3 
4.42115E-3, 4.4128E-3, 
4.45691E-3, 4.48672E-3 
4.51115E-3, 4.5499E-3, 
4.55586E-3, 4.55884E-3 
4.5797E-3, 4.61964E-3, 



1.84488E-3 
2.47911E-3 
2.86418E-3 
3.16818E-3 
3.32376E-3 
3.48649E-3 
3.64683E-3 
3.75174E-3 
3.84414E-3 
3.95084E-3 
4.01819E-3 
4.11595E-3 
4.22086E-3 
4.29835E-3 
4.32994E-3 
4.39194E-3 
4.44678E-3 
4.49268E-3 
4.52904E-3 
4.56361E-3 
4.6107E-3, 
4.68461E-3 
4.72753E-3 
4.78476E-3 
4.79727E-3 
4.83006E-3 
4.87357E-3 
4.89622E-3 
4.90874E-3 
4.93199E-3 
4.97371E-3 
5.02379E-3 
5.0524E-3, 
5.10724E-3 
5.09889E-3 
5.12452E-3 
5.17281E-3 
5.17221E-3 
5.1734E-3, 
5.19725E-3 
5.23063E-3 
5.26639E-3 
5.29321E-3 
5.33434E-3 
5.32242E-3 
.005357, 5 
5.38084E-3 
5.3862E-3, 
5.37965E-3 



4.62262E-3, 4.66494E-3, 
4.67269E-3, 4.68163E-3 



4.70727E-3 
4.76211E-3 
4.81277E-3 
4.8688E-3, 
4.87834E-3 
4.90934E-3 
4.93378E-3 
4.95404E-3 
4.9761E-3, 
5.00352E-3 

5.0369E-3, 
5.07743E-3 
5.11618E-3 
5.16387E-3 
5.16565E-3 
5.17817E-3 

5.20678E-3, 
5.21871E-3 
5.22407E-3 
5.2503E-3, 
5.26937E-3 
5.31169E-3 
5 .34686E-3 

38263E-3, 5 
5.37607E-3 

5.39157E-3, 
5.41243E-3 



4.71203E-3 
4.75078E-3 
4.80383E-3 
4.84973E-3, 
4.89265E-3 
4.94451E-3 
4.94272E-3 
4.98444E-3 



5.01604E-3 

.03392E-3, 
5.0673E-3, 
5.10187E-3 
5.14062E-3 
5.17817E-3 
5.22169E-3 

5.21632E-3, 
5.23838E-3 
5.26222E-3 

5.25745E-3, 
5.27176E-3 
5.30514E-3 
5.32361E-3 

35461E-3, 5 
5.39097E-3 

5.43866E-3, 
5 .42137E-3 



626.7, 

.00302, 6.07389, 

1.40795E-3, 1.68155E-3, 
2.2985E-3, 2.43559E-3, 

2.74794E-3, 2.79682E-3 



3.01201E-3 
3.23256E-3 
3.44953E-3 
3.55087E-3 

.69869E-3, : 
3.82506E-3 
3.94309E-3 
4.04263E-3 
4.10284E-3 
4.1696E-3, 
4 .24292E-3 
4.29716E-3 
4.36571E-3 

4.40505E-3, 
4.46406E-3 

4.51712E-3, 
4.60414E-3 

4.64229E-3, 

4.66316E-3, 
4.72336E-3 
4.74601E-3 
4.76866E-3 
4.82589E-3 

4.82767E-3, 
4.87476E-3 
4.91828E-3 
4.97729E-3 
4.99637E-3 



5.0071E-3, 5.0071E-3, 5 



5.05776E-3 

.0673E-3, 5 

5.08339E-3, 

5.12333E-3 

5.11618E-3 

5.15016E-3 

5.19605E-3 

5.23599E-3, 

5.25685E-3 

5.25924E-3 

5.28368E-3, 

5.31229E-3 

5.30514E-3 

5.3409E-3, 

33375E-3, 5 

5.37667E-3 

5.40766E-3, 

5.43926E-3 



3.07459E-3 
3.31303E-3 
3.45609E-3 
3.58603E-3 

.73088E-3, 
3.82149E-3 
3.94249E-3 
4.03667E-3 
4.12549E-3 

4.19344E-3, 
4.26914E-3 
4.32518E-3 
4.37644E-3 

4.41161E-3, 
4.47062E-3 

4.52248E-3, 
4.58268E-3 

4.62203E-3, 

4.69475E-3, 
4.72813E-3 
4.75555E-3 
4.79847E-3 
4.82469E-3 

4.83125E-3, 
4.88669E-3 
4.90993E-3 
4.9457E-3, 
4.97789E-3 

03988E-3, 
5.05538E-3 

07624E-3, 

5.10903E-3, 
5.1126E-3, 
5.1275E-3, 
5.16983E-3, 
.005174, 

5.20976E-3, 
5.23838E-3, 
5.28725E-3, 

5.29798E-3, 
5.31229E-3, 
5.32838E-3, 

5.33673E-3, 

34984E-3, 
5.38322E-3, 

5.38322E-3, 
5 .41005E-3 
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Appendix III. The Helium Data Reduction Program 

PROGRAM HE (INPUT, OUTPUT, TAPEl) 

C HELIUM COND. DATA REDUCTION, FINAL DR/DT, REVISED 9 MAR 81 
C FOR LS ALONG THE DATA: CALL LS 
C FOR PLOT: CALL PL0T(G(1) ,G(2)) 

DIMENSION DTR(250) ,DTC(250) ,TLOG(250) ,VOF(250) 

DIMENSION F(40),G(10) 

DIMENSION BFY(5) ,GNU(5) ,DT1(250) ,DT2(250) ,DT3(250) ,DT4(250) , 

1 DT5(250),VOLT(250),QI(250) 
DIMENSION ITITLE(13),NAME(8) 
COMMON/DDC/LU,LUC,IFL 

COMMON/8/ IN, IK, P,T,DEN, E,H,S, CV,CP,CSAT, W,WK 

COMMON/LSS/DTR,DTC,TLOG,Q,PRTT 

COMMON /FITT/ F,Y,NFUN 

COMMON /DATATC/R1 ,R2 ,R3Z ,R4Z ,R3B ,R4B , PRESS ,TEMP ,VOF ,TINC 

DATA BFY/0. 509927, -0.338935, 0.271009, -0.232425, 0.206431/ 

DATA GNU/2.4048,5.5201,8.6537,11.7315,14.9309/ 

DATA NAME/10HHANS R0DER,10H X 3528 ,10H , 

110H ,10H ,10H ,10H 

210H / 

C PT, HELIUM, 0.0005 INCH WIRE 

C INITIALIZE 

IFL=156 

CALL GRPH(0,0,8,NAME,0,0) 
Rl=100.1299 
R2=100.1428 
RHOPT=21.37 
A=0. 000635 
B=0.45 
PI=3. 14159 
DENT=20 . 
PRINT 1 
1 F0RMAT(1H1/* RUNPT PRT-T WIRE-T R3Z R4Z TINC 
1VOLTS PRESS DEN T-EXP TC-EXP STAT TC-TN PCT 

2 POWER*/) 

C READ DATATC 

DO 200 K=l,78 
TTRIAL=1.0 

READ(1,*) IDATE,IRUNPT,PRT,PRT1, COUNTS, BAR 
IF(E0F(1)) 9999,100 
100 CONTINUE 

READ(1,*) R3Z,R3B,R4Z,R4B,TINC,VNOVA 

DO 11 1=1,50 

L=(I-1)*5 

READ(1 ,*) V0F(L+1) ,V0F(L+2) ,V0F(L+3) ,V0F(L+4) ,V0F(L+5) 
IF(EOF(l))9999,ll ' urvijT:,; 
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11 CONTINUE 
C 

C INITIAL CALC. AND CORRECTIONS 
C 

DO 12 1=1,250 

12 VOF(I)=VOF(I)-0. 000012 
30 R=PRT/2000. 

PRTT=TEMPP(R) 

TEMP=PRTT 

CPPT=0 . 129+0 .004/80 .*(TEMP-220 . ) 

TCPT=0. 74-0. 03/100. *(TEMP-300.) 

DIFPT=TCPT/RHOPT/CPPT 

PSIA=0. 8728+1. 0351740*COUNTS-0.183212E-6*COUNTS**2 

R3C=-0.319308E-1+0.637332E-3*TEMP+0.151192E-5*TEMP**2 
1 -0.262966E-8*TEMP**3 

R4C=-0 .806515E-1+0 .151734E-2*TEMP+0 .582652E-5*TEMP**2 
1 -0.899257E-8*TEMP**3 

R3Z=R3Z-R3C 

R3B=R3B+R3C 

R4Z=R4Z-R4C 

R4B=R4B+R4C+0 .0024 

PRESS=PSIA/14.696 

VOLT2=2.0*VNOVA 

VEST=0.5 

DO 4 1=1,50 

HPVOLT=VEST+0.1 

VMIN=HPVOLT*0.99 

VMAX=HPVOLT*1.01 

IF(VOLT2.GT.VMIN.AND.VOLT2.LT.VMAX) GO TO 5 

4 VEST=VEST+0.5 

5 CONTINUE 
R3=R3Z+R3B 
R4=R4Z+R4B 

RBRIDGE=1.0/(1.0/(R3+R4)+1.0/(R1+R2)) 
RCIR=100 .3288+RBRIDGE 
VOLTS=HPVOLT/RCIR*RBRIDGE 

C CALCULATE DELTA-T AND T-CORRECTIONS 

IST=51 

IF(IRUNPT.GE.9070) IST=15 

ISTOP=250 

IF(IRUNPT.EQ.9070) ISTOP=220 

IF(IRUNPT.EQ.9071) ISTOP=180 

IF(IRUNPT.EQ.9072) ISTOP=155 

IF(IRUNPT.EQ.9073) ISTOP=126 

IF(IRUNPT.EQ.9074) ISTOP=110 

IF(IRUNPT.EQ.9075) ISTOP= 82 

IF(IRUNPT.EQ.9076) ISTOP= 25 

IF(IRUNPT.EQ.9077) ISTOP= 25 

V250=VNOVA/1 .00319 

V50=V250/1. 00033 

CV2=(V250-V50)/(ALOG(TINC*250)-ALOG(TINC*50)) 

CV1=V50-CV2*AL0G(TINC*50) 
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DO 8 1=1,250 

8 V0LT(I)=CV1+CV2*AL0G(TINC*I) 
DTW=(FINDDT(VOLT(IST) ,V0F(IST) ,TTRIAL)+FINDDT(V0LT(IST0P) , 

1 V0F(IST0P),TTRIAL))/2. 
TEM=DTW+TEMP 
RESW3=FIND1(TEM) 
RESW4=FIND2(TEM) 
RT0TAL=R3B+R4B+RESW3+RESW4 
IQ=(IST+IST0P)/2+l 

Q=(VOLT(IQ)/RTOTAL)**2*(RESW3+RESW4)*100./15.596 
DEN=FNDDHE( PRESS ,TEMP) 
RH0GAS=DEN*4 .0026/1000 . 
TC=THERHE(DEN, TEMP) /10 . 
CPGAS=CPHE( DEN, TEMP) /4 .0026 
CVGAS=CVHE ( DEN , TEMP ) 
DIFGAS=TC/100 ./RHOGAS/CPGAS 
QQ=Q/(4.*PI*TC) 

QQQ=Q*8 .3143*10. /(DEN*CPGAS*39 .948*CVGAS*PI*B**2) 
QQQQ=8.*PI*A/100.*5.6697E-8*TEMP**3/Q 
TIME=0 . 
DO 3 1=1,250 
TIME=TIME+TINC 
TLOG(I)=ALOG(TIME*1000 . ) 
DELT=QQ*AL0G(4 .0*DIFGAS*TIME/A**2/1 .781) 
DTR(I)=FINDDT(VOLT(I) ,V0F(I) ,TTRIAL) 
TTRIAL=DTR(I) 

DTl(I)=A**2*(RH0PT*CPPT-RH0GAS*CPGAS)/(2 .*TC/100 .*TIME)*DELT 

1 -QQ*A**2/(4.*DIFGAS*TIME)*(2.-DIFGAS/DIFPT) 
DT2(I)=0. 

BKT=B**2/(DIFGAS*TIME) 
IF(BKT.GT.5.78) GO TO 10 
BRAK=0. 
DO 9 J=l,5 

9 BRAK=BRAK+EXP(-GNU(J)**2*DIFGAS*TIME/B**2)*(PI*BFY(J))**2 

DT2(I)=QQ*(AL0G(4.0*DIFGAS*TIME/B**2/1.781)+BRAK^ 
10 CONTINUE 

DT3(I)=-QQQ*TIME 
DT4(I)=0. 

DT5(I)=QQQQ*DELT**2 

DTC(I)=DTR(I)+DT1(I)+DT2(I)+DT3(I)+DT4(I)+DT5(I) 
TTT=TEMP+DTR(I) 

R3=FIND1(TTT)+R3B 

R4=FIND2(TTT)+R4B 

RBRIDGE=1.0/(1.0/(R3+R4)+1.0/(R1+R2)) 

VOLT(I)=HPVOLT*RBRIDGE/(RBRIDGE+100.3288) 
RT0TAL=R3+R4 

RESW3=R3-R3B 

RESW4=R4-R4B 

3 00^1^ 
DO 54 1=1,250 
54 DTC(I)=DTC(I)*QI(IQ)/ QI(I) 

LEAST SQUARE DELTA-T VS. LOG(TIMES) 
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: CALL LS 

DO 7 1=1,2 
7 G(I)=0.0 
NFUN=2 
NF=2 

NO PTS=ISTOP-IST+l 
DO 17 L=IST,ISTOP 
F(l)=1.0 
F(2)=TLOG(L) 

16 Y=DTC(L) 

17 CALL FIT 
CALL COEFF 
DO 20 1=1, NF 

20 G(I)=F(I) 

DO 13 1=1,250 

DTR(I)=(DTC(I)-(G(1)+G(2)*TL0G(I)))/DTC(I)*100. 
13 CONTINUE 

CALL PLOTl(IRUNPT) 

CALL STAT 

TCSTAT=F(2)/G(2) 

I CALCULATE THERMAL CONDUCTIVITY 

DEN=FNDDHE ( PRESS ,TEM) 

DENT=DEN 

TC=THERHE(DEN,TEM)/10. 

TCEXP=Q/(4.*3,14159*G(2)) 

PCT=(TCEXP-TC) /TCEXP*100 . 

PRINT 50, IRUNPT,PRTT, TEMP, R3Z,R4Z,TINC, VOLTS, PRESS, DEN, 
1 TEM,TCEXP,TCSTAT,TC,PCT,Q 
50 FORMAT(I6,2F9.3,2F9.4,F9.5,F7.2,F9.3,F8.4,F9.3,F9.5,F8.3, 
1 F9.5,F8.2,F9.5) 

I PLOT RESULTS 

CALL PL0T(G(1),G(2),IRUNPT) 
200 CONTINUE 
9999 CONTINUE 

CALL ENDGRPH(0, 0,0, 0,0,0) 

CALL DDDF 

END 

FUNCTION 0FFSV(EB,RR1,RR2) 
I CALCULATES OFFSET VOLTAGE WHEN RESISTANCES ARE KNOWN 

C0MM0N/DATATC/R1,R2,R3Z,R4Z,R3B,R4B, PRESS, TEMP, VOF,TINC 

DATA(RG=1.0E+07) 

R3=RR1+R3B 

R4=RR2+R4B 

DIV=R1*R2*R3+R2*R3*R4+R3*R4*R1+R4*R1*R2+RG*(R1+R2)*(R3+R4) 

GALVI=EB* (R2*R3-R1*R4 ) /DIV 

OFFSV=GALVI*RG 

RETURN 

END 

FUNCTION FINDl(TIN) 
I WIRE RESISTANCE AS A FUNCTION OF TEMPERATURE 
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C0MM0N/DATATC/R1 ,R2 ,R3Z ,R4Z ,R3B ,R4B , PRESS , TEMP ,VOF ,TINC 

P=PRESS*14.696 

IF(TIN.GT.150.) GO TO 2 

FIND1=-10 .231206+0 .3670809*TIN-0 .9903371E-4*TIN**2-0 .5375358E-5*P 

RETURN 
2 FINDl=-9. 0654718+0. 35344447*TIN-0.59234427E-4*TIN**2-0.966271E-5*P 

RETURN 

END 

FUNCTION FIND2(TIN) 

WIRE RESISTANCE AS A FUNCTION OF TEMPERATURE 

COMMON/DATATC/Rl^jRSZ^ZjRSB^B, PRESS, TEMP, VOF,TINC 

P=PRESS*14.696 

IF(TIN.GT.150.) GO TO 2 

FIND2=-5 .057558+0 .1823527*TIN-0 .525728E-4*TIN**2-0 .285743E-5*P 

RETURN 
2 FIND2=-4 .346459+0 .17402506*TIN-0 .2831553E-4*TIN**2-0 ,452695E-5*P 

RETURN 

END 

FUNCTION FINDDT(V1,V2,T1) 

AN APPROXIMATION, TO REPLACE WIRE-R 

COMMON /DATATC/R1 ,R2 ,R3Z ,R4Z ,R3B ,R4B ,PRESS , TEMP , VOF ,T INC 

V=V2 

TT=T1+TEMP 

DO 10 1=1,10 

RE1=FIND1(TT) 

RE2=FIND2(TT) 

V0FF1 =0FFSV(V1,RE1,RE2) 

IF(ABS(V-V0FFl)-1.0E-6 ) 20,20,1 
1 T=TT+0.01 

RE1=FIND1(T ) 

RE2=FIND2(T ) 

VOFF2 =0FFSV(V1,RE1,RE2) 

DVDT=(VOFF2-VOFFl) /0 .01 

CORR= (V-V0FF1 ) /DVDT 
10 TT=TT+CORR 
20 FINDDT=TT-TEMP 

RETURN 

END 

FUNCTION TEMPP(R) 

ABBR. FOR 1776265 FROM 70 TO 310 K, CORR. 12/19/80 

DIMENSION RR(241) 

4.10575, 4.21531, 4.32508, 4.43505, 
4.87649, 4.98714, 5.09788, 5.20868, 
5.65235, 5.76332, 5.87430, 5.98527, 
6.31807, 6.42895, 6.53979, 6.65058, 
7.09327, 7.20379, 

7.42465, 7.53498, 7.64523, 7.75541, 

8.19538, 8.30517, 8.41489, 8.52453, 

8.96226, 9.07149, 9.18065, 9.28971, 

9.72519, 9.83386, 9.94245,10.05096, 

,10. 48424, 10. 5923/, 10. 70043, 10. 80841, 

11.23962,11.34725,11.45480,11.56229, 

11.99158,12.09873,12.20583,12.31286, 

,12 .74036 ,12 .84708 , 12 . 95375 ,13 .06036 , 



DATA(RR(I) ,1=1 ,131) /3 .99645 

4.54520, 4.65550, 4.76594 

5.31954, 5.43045, 5.54139 

6.09622, 6.20716 

6.76134, 6.87204, 6.98268 

7.31426 

7.86552, 7.97555, 8.08550 

8.63408, 8.74356, 8.85295 

9.39870, 9.50761, 9.61644 

10 .15939 ,10 .26775 ,10 .37603 

10.91632,11.02416,11.13192 

11.66971,11.77707,11.88435 

12.41982,12.52673,12.63357 
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H 13.16690,13.27340,13.37983,13.48621,13.59254,13.69881,13.80502, 
I 13.91118,14.01729,14.12335,14.22935,14.33530,14.44121,14.54706, 
J 14.65286,14.75861,14.86431,14.96997,15.07557,15.18113,15.28664, 
K 15.39210,15.49752,15.60289,15.70822,15.81350,15.91873,16.02392, 
L 16.12907,16.23417,16.33923,16.44424,16.54922,16.65414,16.75903, 
M 16.86387,16.96868,17.07344,17.17816,17.28284,17.38747,17.49207, 
N 17.59663,17.70114,17.80562,17.91006,18.01445/ 

DATA(RR(I), 1=132, 241)/ 18.11881,18.22313,18.32741,18.43165, 
A 18.53586,18.64002,18.74415,18.84824,18.95229,19.05630,19.16028, 
B 19.26422,19.36812,19.47199,19.57582,19.67961,19.78337,19.88709, 
C 19.99078,20.09442,20.19804,20.30162,20.40516,20.50867,20.61214, 
D 20.71558,20.81898,20.92235,21.02568,21.12898,21.23225,21.33548, 
E 21.43867,21.54184,21.64497,21.74806,21.85112,21.95415,22.05715, 
F 22.16011,22.26304,22.36594,22.46880,22.57163,22.67443,22.77720, 
G 22 .87993 ,22 .98263 ,23 .08530,23 .18793 ,23 .29054 ,23 .39311 ,23 .49565 , 
H 23 .59816 , 23 . 70063 ,23 .80308 ,23 . 90549 ,24 .00787 ,24 .11022 , 24 .21254 , 
I 24 .31483,24 .41708,24 .51931 ,24 .62150 ,24 .72366 ,24 .82579 ,24 .92789 , 
J 25.02996,25.13199,25.23400,25.33597,25.43791,25.53982,25.64170, 
K 25.74355,25.84537,25.94715,26.04891,26.15063,26.25232,26.35398, 
L 26.45561,26.55721,26.65877,26.76031,26.86181,26.96329,27.06473, 
M 27.16614 ,27.26752,27.36887,27.47019,27.57147,27.67273, 

1 27.77395,27.87514,27.97631,28.07744,28.17854,28.27961,28.38065, 

2 28 .48165 ,28 .58263 ,28 .68357 ,28 .78449 ,28 .88537 ,28 .98622 ,29 .08704 , 

3 29.18783,29.28859/ 
DO 3 1=1,241 
IF(R.LT.RR(I)) GO TO 4 

3 CONTINUE 

4 TEMPP=I 
TEMPP=TEMPP+68.+(R-RR(I-l))/(RR(I)-RR(I-l)) 
RETURN 

END 
SUBROUTINE LS 

LEAST SQUARE ALONG THE DATA 

DIMENSION DTR(250) ,DTC(250) ,TL0G(250) ,VOF(250) 

DIMENSION F(40),G(10) 

DIMENSION RES(20) ,ERR(20) ,IM(20) ,TM(20) 

COMMON/LSS/DTR,DTC,TLOG,Q,PRTT 

COMMON /FITT/ F,Y,NFUN 

IST=11 

DO 13 J=l,16 

ISTOP=IST+89 

IM(J)=(IST+IST0P)/2 

TM(J)=PRTT+DTC(IM(J)) 

NFUN=2 

NF=2 

NO PTS=IST0P-IST+1 

DO 17 L=IST,ISTOP 

F(l)=1.0 

F(2)=TL0G(L) 

16 Y=DTC(L) 

17 CALL FIT 
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CALL COEFF 

RES(J)=Q/(4.*3.14159*F(2)) 
CALL STAT 
ERR(J)=RES(J)/(Q/(4.*3.14159*F(2))) 

13 IST=IST+10 

PRINT 15,((IM(J),TM(J)),J=1,10) 
PRINT 14,((RES(J),ERR(J)),J=1,10) 
PRINT 15,((IM(J),TM(J)),J=11,16) 
PRINT 14,((RES(J),ERR(J)),J=11,16) 

14 FORMAT (5X,20F6. 4) 

15 F0RMAT(5X,10(I4,F8.3)) 
RETURN 

END 

SUBROUTINE PLOT (Gl , G2 , IDNO ) 

PLOTS RESULTS 

DIMENSION ITITLE(13) ,NAME(8) ,ICOMM(2) 

DIMENSION DTR(250) ,DTC(250) ,TLOG(250) ,VOF(250) 

COMMON/LSS/DTR,DTC,TLOG,Q,PRTT 

DATA ITITLE/10H ,10H HELIUM T,10HHERMAL CON, 
110HDUCTIVITY ,10H ,10H ,10HLN(T) T IN, 

21 OH MILLISEC ,10H ,10H ,10HDELTA T, , 

310HT IN K ,10H / 

DATA ICOMM(l)/10H RUNPT / 

ENCODE(10,l,ICOMM(2))IDNO 
1 FORMAT(IIO) 

NO PTS=240 

DO 51 1=11,250 

L=I-10 

TLOG(L)=TLOG(I) 
51 DTC(L)=DTC(I) 

CALL PGRPH(TLOG,DTC,NO PTS,ITITLE,0,6) 

DTC(1)=G1+G2*TL0G(1) 

DTC(2)=G1+G2*TL0G(N0 PTS) 

TL0G(2)=TL0G(N0 PTS) 

CALL CLGRPH(TL0G,DTC,2,0,0,0) 

ABFRACT=0.05 

ORDFR=0.9 

CALL COMGRPH( ABFRACT , ORDFR , 2 , ICOMM) 

RETURN 

END 

SUBROUTINE PLOTl(IDNO) 

PLOTS DEV. FROM LS LINES 

DIMENSION ITITLE(13) ,NAME(8) ,TT(250) ,IC0MM(2) 

DIMENSION DTR(250) ,DTC(250) ,TLOG(250) ,VOF(250) 

COMMON/LSS/DTR,DTC,TLOG,Q,PRTT 

DATA ICOMM(1)/10H RUNPT / 

ENC0DE(10,1,IC0MM(2))IDN0 
1 FORMAT(IIO) 

NO PTS=240 

DO 51 1=11,250 
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L=I-10 

TT(L)=TLOG(I) 
51 DTR(L)=DTR(I) 

CALL PGRPH(TT,DTR,NO PTS,ITITLE,0,6) 

ABFRACT=0.05 

0RDFR=0.9 

CALL C0MGRPH(ABFRACT,0RDFR,2,IC0MM) 

RETURN 

END 

SUBROUTINE FITTER 
C CHANGED 13 SEPT., 1973 OLDER VERSION OF FITTER 

C COMMON STATEMENT ADDED TO MATCH CALLING SEQUENCE OF NEWER VERSION 
C SPECIAL, TO GET RESS IN F(10) 

DIMENSION F(40),A(40,41) 

COMMON /FITT/ F,Y,NFUN 

DOUBLE PRECISION A,SY,SYY,RES,DET 

EQUIVALENCE ( NC , FNC ) 

DATA (NTR=-1) $ DATA(NDIM=40) 

ENTRY FIT 

IF(NTR) 1,3,3 

1 NP=0 
NF=NFUN 

IF(NF.GT.NDIM) GO TO 44 
NC=0 

SY=0. 
SYY=0 . 
NY=NF+1 
DO 2 1=1, NY 
DO 2 J=1,NF 

2 A(J,I)=0. 
IF(NTR.EQ.O) GO TO 11 
NTR=0 

3 SY=SY+Y 
SYY=SYY+Y*Y 
DO 4 J=1,NF 
A(J,NY)=A(J,NY)+Y*F(J) 
DO 4 1=1, NF 

4 A(I,J)=A(I,J)+F(I)*F(J) 
NP=NP+1 

RETURN 
ENTRY CONSTR 
IF(NTR) 10,11,11 

10 NTR=0 
GO TO 1 

11 N=NY-1 
IF(NY+NC+2.GT.NDIM) GO TO 44 
DO 12 1=1 ,N 
A(I,NY+1)=A(I,NY) 

12 A(I,NY)=A(NY,I)=A(NDIM-NC,I)=F(I) 
NC=NC+1 

DO 13 I=NF,N 

13 A(I+1,NY)=A(NY,I+1)=0. 
NY=NY+1 
A(NY-1,NY)=Y 
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RETURN 

ENTRY COEFF 

N=NY-1 

DO 20 1=1, NF 

20 F(I)=A(I,NY) 
IF(N.EQ.1)G0 TO 23 
DO 22 1=2, N 

DO 21 J=I,NY 

21 A(I-1,J)=A(I-1,J)/A(I-1,I-1) 
DO 22 J=I,N 

DO 22 K=I,NY 

22 A(J,K)=A(J,K)-A(J,I-1)*A(I-1,K) 

23 A(N,NY)=A(N,NY)/A(N,N) 
IF(N.EQ.1)G0 TO 26 

DO 24 1=2, N 

L=N-I+2 

DO 24 J=L,N 

24 A(L-1 ,NY)=A(L-1 ,NY)-A(L-1 , J) *A( J ,NY) 

26 NTR=-1 
RES=SYY 

DO 25 1=1, NF 
RES=RES-F(I)*A(I,NY) 

25 F(I)=A(I,NY) 
NFUN=NP 
DF=NP-NF+NC 
Y=FNC 
FNP=NP 
RETURN 
ENTRY STAT 
IF(NC.LE.O)GO TO 28 
DO 27 1=1, NC 

DO 27 J=1,NF 

27 RES=RES-A(NDIM-I+1,J)*A(J,NY)*A(NF+I,NY) 

28 TOT=SYY-SY*SY/FNP 
REG=TOT-RES 
SYY=RES/DF 

ST=1.96+2.72/DF+8.04/DF**3 
DET=1.0D+0 
DO 30 1=1, NF 
DET=DET*A(I,I) 

30 A(I,I)=1.0D+0/A(I,I) 
IF(NF.EQ.1)G0 TO 322 
DO 32 1=2, NF 

DO 32 J=2,I 

SY=0. 

DO 31 K=J,I 

31 SY=SY-A(I,K-1)*A(K-1,J-1) 

32 A(I,J-1)= S Y*A(I,I) 
322 CONTINUE 

DO 36 1=1 ,NF 

L=NF-I 

DO 33 J=1,L 

K=NF-J 

DO 33 M=1,J 

492 



N=NF-M+1 

33 A(K,I)=A(K,I)-A(K,N)*A(N,I) 
IF(I.EQ.1)G0 TO 344 

DO 34 J=2,I 

34 A(J-1,I)=A(I,J-1)*SYY 
344 CONTINUE 

DO 35 J=1,I 

35 A(I,J)=A(I,J)*SYY 
IF(A(I,I).LE.O.O)A(I,I)=-A(I,I) 

36 F(I)=ST*DSQRT(A(I,I)) 

PRINT 37,(A(I,NY),F(I),I=1,NF) 

37 FORMAT(*lTHE COEFFICIENTS AND THEIR ESTIMATED ERRORS ARE:*/ 
1(* *E18.10,* +OR-*E9.2)) 

SYY=DSQRT(SYY) 

PRINT 38,RES,REG,TOT,SYY,DET 

F(10)=RES 

38 FORMAT( * ESTIMATED RESIDUAL SUM OF SQUARES =*E17.9/ 

1 * ESTIMATED REGRESSION SUM OF SQUARES =*E17.9/ 

2 * ESTIMATED TOTAL SUM OF SQUARES =*E17.9/ 

3* STANDARD DEVIATION =*E17.9/* DETERMINANT OF THE MATRIX =*E17.9/ 
4 * DISPERSION MATRIX* ) 
DO 39 J=1,NF 

39 PRINT 40,(A(I,J),I=1,NF) 

39 CONTINUE 

40 FORMAT(7E19.9) 
PRINT 46 

46 F0RMAT(1H ) 
RESSS=RES 
Y=SQRT(RESSS/DF) 
NFUN=DF 

43 RETURN 

44 PRINT 45 

45 F0RMAT(*1THE ARRAYS IN THE FITTING PROGRAM ARE TOO SMALL TO HOLD T 
1HE NUMBER OF CONSTRAINTS AND FUNCTIONS ASKED FOR IN THE CALLING PR 
20GRAM*) 

STOP 
END 
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